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A low momentum nucleon-nucleon (NN) potential Vi ow -k is derived from meson exhange po- 
tentials by integrating out the model dependent high momentum modes of Vnn- The smooth 
and approximately unique Vi ow -k is used as input for shell model calculations instead of the usual 

Q-t. 



Brueckner G matrix. Such an approach eliminates the nuclear mass dependence of the input interac- 
tion one finds in the G matrix approach, allowing the same input interaction to be used in different 
. nuclear regions. Shell model calculations of 18 0, 134 Te and 135 I using the same input Vi ow ~k have 

^vq been performed. For cut-off momentum A in the vicinity of 2 /m _1 , our calculated low-lying spectra 

for these nuclei are in good agreement with experiments, and are weakly dependent on A. 
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■ 21.60.Cs; 21.30.Fe; 27.80.+j 
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, A fundamental problem in nuclear physics has been the determination of the effective nucleon-nucleon (NN) inter- 
action used in the nuclear shell model, which has been successful in describing a variety of nuclear properties. There 
have been a number of successful approaches [1-4] for this determination, ranging from empirical fits of experimental 
data, to deriving it microscopically from the bare NN potential. Despite impressive quantitative successes, the tradi- 
tional microscopic approach suffers the fate of being " model dependent" owing to the fact that there is no unique Vnn 
to start from. Moreover, as the Brueckner G matrix has traditionally been the starting point, one obtains different 
input interactions for nuclei in different mass regions as a result of the Pauli blocking operator. 

In this work, we propose a different approach to shell model effective interactions that is motivated by the recent 
' applications of effective field theory (EFT) and the renormalization group (RG) to low energy nuclear systems [5-8] . 
Our aim is to remove some of the model dependence that arises at short distances in the various Vnn models, and also 
to eliminate the mass dependence one finds in the G matrix approach, thus allowing the same interaction to be used 
in different nuclear regions such as those for 18 and 134 Te. A central theme of the RG-EFT approach is that physics 
in the infrared region is insensitive to the details of the short distance dynamics. One can therefore have infinitely 
many theories that differ substantially at small distances, but still give the same low energy physics if they possess 
the same symmetries and the "correct" long- wavelength structure [5,8]. The fact that the various meson models for 
Vnn share the same one pion tail, but differ significantly in how they treat the shorter distance pieces illustrates this 
explicitly as they give the same phase shifts and deuteron binding energy. In RG language, the short distance pieces 
of Vn n are like irrelevant operators since their detailed form can not be resolved from low energy data. 

Motivated by these observations, we would like to derive a low-momentum NN potential Vi ow -k by integrating 
out the high momentum components of different models of Vnn m the sense of the RG [5,8], and investigate its 
suitability of being used directly as a model independent effective interaction for shell model calculations. We shall 
use in the present work the CD-Bonn NN potential [9] for Vnn- In the following, we shall first describe our method 
for carrying out the high-momentum integration. Shell model calculations for 18 0, 134 Te and 135 I using Vi ow -k will 
then be performed. Our results will be discussed, especially about their dependence on the cut-off momentum A. 

The first step in our approach is to integrate out the model dependent high momentum components of Vnn- In 
accordance with the general definition of a renormalization group transformation, the decimation must be such that 
low energy observables calculated in the full theory are exactly preserved by the effective theory. We turn to the 
model space methods of nuclear structure theory for guidance, as there has been much work in recent years discussing 
their similarity to the Wilson RG approach [4,10,11]. While the technical details differ, both approaches attempt 
to thin-out, or limit the degrees of freedom one must explicitly consider to describe the physics in some low energy 
regime. Once the relevant low energy modes are identified, all remaining modes or states are "integrated" out. Their 
effects are then implicitly buried inside the effective interaction in a manner that leaves the low energy observables 
invariant. One successful model-space reduction method is the Kuo-Lee-Ratcliff (KLR) folded diagram theory [12,13]. 
For the nucleon-nucleon problem in vacuum, the RG approach simply means that the low momentum T matrix and 
the deuteron binding energy calculated from Vnn must be reproduced by Vi ow -k, but with all loop integrals cut off 
at some A. Therefore, we start from the half-on-shell T-matrix 

f°° 1 

T(k',k,k 2 ) = V NN (k' 1 k) + q 2 dqV N N(k',q)- 2 — — -T(q, k, k 2 ). (1) 
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We then define an effective tow-momentum T-matrix by 



f A 1 
Ti ow -k{p',P,p 2 ) = Vi ow - k (p',p) + / q dqVi ow -k{p',q) 2 _ 2 .„ + Ti ow -k{q,P,P ), (2) 

where A denotes a momentum space cut-off (such as A=2/to~ 1 ) and (p',p) < A. We require the above T-matrices 
satisfying the condition 

T(p',p,p 2 ) = T low - k (p\p,p 2 ); (p\p) < A. (3) 

The above equations define the effective low momentum interaction Viow-k- In the following, let us show that the 
above equations are satisfied by the solution 

V low - k = Q-Q> J Q + Q> J qJq-Q'JqJqJq+... , (4) 

which is just the KLR folded-diagram effective interaction [12,13]. A preliminary account of this result has been 
reported as a work in progress at a recent conference [14]. 

In time dependent formulation, the T-matrix of Eq.(l) can be written as {k' \ VU(0, — oo) k), U being the time 
evolution operator. In this way we can readily perform a diagrammatic analysis of the T-matrix. A general term of 
it may be written as {k' \ (V + V-^^V + V-^V-^^V + ■ ■ •) | k) where e(k) = (k 2 — H Q ), H Q being the unperrurbed 

Hamiltonian. Note that the intermediate states (represented by 1 in the numerator) cover the entire space, and 
1 = P + Q where P denotes the model space (momentum < A) and Q its complement. Expanding it out in terms of P 
and Q, a typical term of T is of the form V^V^V^V^V^V. Let us define a Q-box as Q = V+V^V+V^V§V+-, 
where all intermediate states belong to Q. One readily sees that the T-matrix can be regrouped as a Q-box series, 
namely (p' \ T \ p) = (p' \ [Q + Q^Q + Q^-Q^-Q + ■ ■ ■] \ p). Note that all the Q-boxes have the same energy variable, 
namely p 2 . 

This regrouping is depicted in Fig. 1, where each Q-box is denoted by a circle and the solid line represents the 
propagator The diagrams A, B and C are respectively the one- and two- and three-Q-box terms of T, and clearly 
T=A+B+C+- • •. Note the dashed vertical line is not a propagator; it is just a "ghost" line to indicate the external 
indices. We now perform a folded-diagram factorization for the T-matrix, following closely the KLR folded-diagram 
method [12,13]. Diagram B of Fig. 1 is factorized into the product of two parts (see Bl) where the time integrations 
of the two parts are independent from each other, each integrating from — oo to 0. In this way we have introduced 
a time- incorrect contribution which must be corrected. In other words B is not equal to Bl, rather it is equal to Bl 
plus the folded-diagram correction B2. Note that the integral sign represents a generalized folding [12,13]. 

Similarly we factorize the three-Q-box term C as shown in the third line of Fig. 1. Higher-order Q-box terms are 
also factorized following the same folded-diagram procedure. Let us now collecting terms in the figure in a "slanted" 
way. The sum of terms Al, B2, C3... is just the low-momentum effective interaction of Eq.(4). (Note that the 
leading Q-box of any folded term must be at least second order in Vnn, and hence it is denoted as Q'-box which 
equals to Q-box with terms first-order in Vnn subtracted.) The sum Bl, C2, D3.... is Vi ow - k ^Q. Similarly the 

sum C1+D2+E3-I is just Vi ow -k^Q^Q- (Note diagrams Dl, D2, • • •, El, E2, • • • are not shown in the figure.) 

Continuing this way, it is easy to see that Eqs. (1) to (3) arc satisfied by the low momentum effective interaction of 
Eq.(4). 

The effective interaction of Eq.(4) can be calculated using iteration methods. A number of such iteration methods 
have been developed; the Krenciglowa-Kuo [15] and the Lee-Suzuki iteration methods [16] are two examples. These 
methods were formulated primarily for the case of degenerate PH P, H being the unperturbed Hamiltonian. For 
our present two-nucleon problem, PH n P is obviously non-degenerate. Non-degenerate iteration methods [17] are 
more complicated. However, a recent iteration method developed by Andreozzi [18] is particularly efficient for the 
non-degenerate case. This method shall be referred to as the Andreozzi-Lee- Suzuki (ALS) iteration method, and has 
been employed in the present work. 

We have carried out numerical checks to ensure that certain low-energy physics of Vnn are indeed preserved by 
Vi ow -k- We first check the deuteron binding energy BEd- We have calculated BEd using Vi ow -k for many values 
of A, and for all cases the BEd given by Vi ow -k agrees very accurately (to 4 places after the decimal) with that 
given by Vnn- (Note that when A approaches oo V ow -k is the same as Vnn-) In Fig. 2, we present some 1 So and 
3 Po phase shifts calculated from the CD-Bonn Vnn (dotted line) and the V ow -k (circles) derived from it, using a 
momentum cut-off A = 2.0/m -1 . As seen, the phase shifts from the former are well reproduced by the latter. We 
have also checked the half-on-shell T-matrix given by Vnn and by Vi ow ~k, and found very good agreement between 
them [14]. The above agreements are expected, as we have shown that the T-matrix equivalence of Eq.(3) holds for 



any A. In short, our numerical checks have reaffirmed that the deuteron binding energy, low energy phase shifts and 
low momentum half-on-shell T-matrix of Vnn are all preserved by Vi ow -k- As far as those physical quantities are 
concerned, Vi ow -k and Vnn are equivalent. 

Having proven the "physical equivalence" of Vi ow -k and Vnn in the sense of the RG, we turn now to microscopic 
shell model calculations in which we use Vi ow -k as the input interaction. A folded-diagram formulation [13,2,3] is 
employed. An important feature here is that this formalism allows us to calculate the energy differences of neighboring 
many body systems. For example, we can calculate the energy difference of ls O and the ground state energy of le O, 
starting from the experimental 17 O single particle energies (s.p.e.) and a shell model effective interaction V e ff derived 
microscopically from an underlying NN potential. This folded diagram method has been rather successfully applied 
to many nuclei using G-matrix interactions. [2,3] There the basic input to the calculation are the matrix elements 
(ni7i2 | G | 113114) where G is the Brueckner G-matrix and the n's are harmonic oscillator wave functions. 

Since Vi ow -k is already a smooth potential, it is no longer necessary to first calculate the G-matrix. Thus in 
our present work, the starting basic input are just the matrix elements {111122 | Viow-k \ 113114), and thereafter 
our calculation procedures are exactly the same as described in references [2,3]. A model space with two va- 
lence neutrons in the (0<i5/2, 0<i3/2, ls^) shell is used for 18 0, and one with two and three valence protons in the 
(0(^7/2, lc?5/2, 1^3/2, 2si/ 2 , 0/in/ 2 ) shell for 134 Te and 135 I, respectively. As customary, we use s.p.e. extracted from the 
experimental spectra of the corresponding single-particle valence nuclei, 17 and 133 Sb [19]. For the absolute scaling 
of the sets of s.p.e., the mass-excess values for 17 and 133 Sb have been taken from Ref. [20,21]. For 134 Te and 135 I, 
we assume that the contribution of the Coulomb interaction between valence protons is equal to the matrix element 
of the Coulomb force between the states (#z) % =0 +- 

As shown in Fig. 3, our calculated low-lying J T states of 18 agree highly satisfactorily with experiments [19]. In 
the same figure, results of the corresponding G-matrix calculations are also shown; the Vi ow -k results are just as good 
or slightly better. It may be mentioned that our Vi ow -k is slightly non-hermitian. A hermitian Vi ow ~k can be obtained 
using the Okubo transformation [7]. We have constructed such a hermitian Vi ow -k using the Suzuki-Okamoto method 
[22]. We have found that the shell model energy levels given by the two Vi ow -kS are very similar, probably because 
our Vi ow -k is only slightly non-hermitian. In a concurrent paper [11] we have found that Vi ow -k is almost independent 
of the underlying Vnn for the values of A considered here. Therefore, although the CD-Bonn potential [9] is used in 
our present calculations, we stress that similar results will be obtained if we calculate Vi ow -k from other models such 
as the Paris or Argonne V-18 potentials. 

It may be mentioned that the G-matrix is energy dependent and Pauli blocking dependent, while Vi ow -k is not. 
This is a desirable feature, indicating that Vi ow ~k may be suitable also for other nuclear regions. To study this point, 
we have used the same Vi ow -k in a shell model calculation of 134 Te as mentioned earlier. It is encouraging that our 
results for 134 Te also agree well with experiments [19] as shown in Fig. 4. Again the Vi ow -k results are just as good 
or slightly better than the G-matrix results. We emphasize that we have used the same Vi ow -k interaction in both 
18 and 134 Te calculations, and it appears to work equally well for both nuclei. This is in marked contrast to the 
traditional approach in which one has to use different G-matrices for different mass regions, as the associated Pauli 
blocking operators are different. This is an appealing result, as it suggests the possibility for a common shell- model 
interaction that attenuates much of the dependence on the Vnn model and is suitable for a wide range of nuclei. 

It is of interest to investigate if the same Vi ow -u is suitable for nuclei with more than two valence nucleons. It is 
primarily for this purpose we have carried out the shell model calculation of 135 I mentioned earlier, using the same 
Vi ow -k- Our results are shown in Fig. 5. It is gratifying that the calculated excitation spectra are in very good 
agreement with experiments [19]. We note that the valence interaction energy for the three valence nucleons given by 
our calculation is slightly overbound, by about 0.3 MeV. This may be an indication of the need of a weak three-body 
force for this nucleus, which has three valence nucleons . (Our Vi ow ~k is a two-body interaction.) We plan to study 
this topic in a future work. 

An important issue is what value one should use for A. Guided by general EFT arguments, the minimum value for 
A must be large enough so that Vi ow ^u explicitly contains the necessary degrees of freedom for the physical system. 
For example, 2tt exchanges are important for low energy nuclear physics, and to adequately include the corresponding 
degree of freedom we need to have A m i n larger than ~ ni2^, i.e. ~ 1.4/m . In fact we have found that Vi ow -k varies 
strongly with A when it is smaller than that value. A general signal for A m i n is when the calculated physical quantities 
first become insensitive to A [5] . Conversely, we want A to be smaller than the short distance scale A max at which the 
model dependence of the different Vnn starts to creep in [5]. Systems in which these two constraints are consistent 
with each other (i.e., A min < A max ) are amenable to EFT-RG inspired effective theories, as they possess a clear 
separation of scales between the relevant long wavelength modes and the model dependent short distance structure. 
We have found [11] that A max should not be much greater than 2.0-2.5 fin -1 as this is the scale after which Vi ow -k 
first becomes strongly dependent on the particular Vnn used. There is another consideration: Most NN potentials 
are constructed to fit empirical phase shifts up to Ei a b « 350 MeV [9]. Since Ei a b < 2h 2 A 2 /M, M being the nucleon 



mass, and one should require Vi ow -k to reproduce the same empirical phase shifts, a choice of A in the vicinity of 2 
fm^ 1 would seem to be appropriate . 

Guided by the above considerations, we have used in our calculations two values for the momentum cut-off, namely 
A = 2.0 and 2.2 fm _1 as shown in Figs. 3 to 5. It is satisfying to see that the results are rather insensitive to the choice 
of A, in harmony with the EFT philosophy mentioned earlier. Perhaps more importantly, both are in satisfactory 
agreement with experiments. 

In summary, we have investigated a RG-EFT inspired approach to shell model calculations that is a "first step" 
towards a model independent calculation that uses one common interaction over a wide range of nuclei. Using the KLR 
folded diagram approach in conjunction with the ALS iteration method, we have performed a RG decimation where 
the model dependent pieces of Vnn models are integrated out to obtain a nearly unique low momentum potential 
Viow-k- This Viow-k preserves the deuteron pole as well as the low energy phase shifts and half-on-shell T matrix. 
We have used VJ ou) _fc, which is a smooth potential, directly in shell model calculations of 18 O, 134 Te and 135 I without 
first calculating the G matrix. The results are all in satisfactory agreement with experiment, and they are insensitive 
to A in the neighborhood of A w 2 /to -1 . We do feel that Vi ow -k may become a promising and reliable effective 
interaction for shell model calculations of few valence nucleons, over a wide range of nuclear regions. 
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FIG 1. Folded-diagram factorization of the half-on-shell T-matrix. 



FIG. 2. Comparison of phase shifts given by Vi ow -k and Vnn- 



FIG. 3. Low-lying states of 18 O. 
FIG. 4. Low-lying states of 134 Te. 
FIG. 5. Low-lying states of 135 I. 
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